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Abstract 

The propagation of elastic waves inside the Earth provides us with in- 
formations about the geological structure of the Earth's interior. Since the 
beginning of seismology, people have been using waves created by earthquakes 
or by artificial explosions. They record the waves as functions of time using 
seismometers located at different stations on the Earth's surface. Even with- 
out any earthquake or explosion, a weak signal is still recorded which has no 
evident structure: it is a "noise". How to use these noises? This is the goal 
of the method of "passive imaging". The main observation is the following 
one: the time correlation of the noisy fields, computed from the fields recorded 
at the points A and B, is "close" to the Green's function G(r, A, B) of the 
wave propagation. The aim of this paper is to provide a mathematical context 
for this approach and to show, in particular, how the methods of semiclassi- 
cal analysis can be be used in order to find the asymptotic behaviour of the 
correlations. 

Introduction 

Seismologists want to recover the physical parameters of the Earth's interior from 
records (called scismograms) , at the Earth's surface, of the elastic waves propagating 
in its interior. From the mathematical point of view, it is an example of a so-called 
"inverse problem": the most famous inverse problems are the Calderon problem 
(recovering the conductance of a domain from boundary measurements) and the 
Kac problem (recovering a Euclidean domain from the spectrum of its Dirichlct 
Laplacian). The use of elastic waves created by an earthquake (resp. an artificial 
explosion) has been quite successful in order to recover the large scale structure 
(resp. the small scale structure, e.g. for oil detection). The geological structure of 
the Earth crust up to depths of 30-50 kilometres is much more difficult to know. 

On the other hand, only a small part of the scismograms was used: this part 
corresponds to the propagation of well identified body and surface waves: typically 
the P-waves, the S- waves and the Rayleigh waves. The last part (called the "coda") 
of the seismograms were not used. This was a long standing question asked by 
K. Aki, one of the founders of modern seismology. Quite recently, it was observed 
by M. Campillo and A. Paul [5] that the time-correlation of coda waves is closely 
related to the Green's function. 

Later, M. Campillo and his collaborators showed that the correlations of ambient 
noise already yield the Green's function [5J [231 [55] ■ Similar facts had been known 
before in acoustics [TFJ [3TJ [351 [5] and helio-seismology [15] . Several theoretical 
results were derived in [51 [5T1 [53J . Applications to clock synchronisation [55] and to 
volcanic eruptions forecasting [4] have been developed. The paper [16] is a recent 
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survey. So, let us assume that we can recover the Green's function: how does it 
help to solve our inverse problem? Here enters scmiclassics (or ray theory): the 
scmiclassical behaviour of the Green's function contains the ray dynamics. Even 
more, the full Green's function is not really needed for that, its phase is sufficient. So 
it is natural to study the problem of passive imaging using the tools of scmiclassical 
analysis. 

The goal of my paper is to present some mathematical models for this method 
of passive imaging. They will be applicable to any kind of wave, but I have in 
mind seismic waves. Semiclassical analysis will be used as a way to get a geometric 
approximation to the wave propagation AND in the modelling of the source noise. 
I do not at all pretend that the assumptions of my models are realistic from the 
physical point of view. I just hope that they will help seismologists by providing 
another point of view and other words. On the other hand, these models put light 
onto several interesting mathematical problems: mode conversion, inverse spectral 
problems, mode-ling of random fields . . . 

Things work as follows: let us assume that we have a medium X (a smooth 
manifold) and a smooth, deterministic (no randomness in it) linear wave equation 
in X. We hope to recover (part of) the geometry of X from the wave propagation 
on it. We assume that there is somewhere in X a source of noise f (x,t) which is a 
stationary random field. This source generates, by the (linear) wave propagation, a 
field u(x,t) (assumed to be scalar in this introduction). This field u is recorded at 
different points A, B, ■ ■ ■ on long time intervals. We want to get some information 
on the propagation of waves from B to A in X from the correlation 



which can be computed numerically from the fields recorded at A and B. In seis- 
mology the medium X will be the Earth, the waves are elastic waves, the main 
source noise is due to the interaction between the ocean or the atmosphere with the 
Earth's crust, and the field is recorded at the Earth's surface. 

It turns out that Ca,b (t) is closely related to the deterministic Green's junction 
of the wave equation in X. It means that one can hope to recover, using Fourier 
analysis, the propagation speeds of waves between A and B as a function of the 
frequency, and, more precisely, the Hamiltonian or, more generally, the so-called 
dispersion relation. 

If the wave dynamics is time reversal symmetric, the correlation admits also a 
symmetry by changing r into — r; this observation has been used for clock synchro- 
nisation, see [22] . 

The goal of this paper is to give precise formulae for Ca,b{ t ) hi the high fre- 
quency limit, under some assumptions on the source f. We will assume that the 
source field is a stationary ergodic random field, the covariance of which K(x, y,s — 
s') = E(/(x, s)f(y, s')) is given by the Schwartz kernel of a pseudodiffercntial oper- 
ator. This is a strong assumption, which implies in particular a rapid decay of the 
correlations of the source f away from the diagonal. 

More precisely, we have two small parameters, one of them entering into the 
decorrelation distance of the source noise, the other one in the high frequency prop- 
agation (the ray method). The fact that both are of the same order of magnitude 
is crucial for the method. As a result, the method works well in certain frequency 
ranges of the wave propagation. From the physical point of view, the main result 
is that, with some assumptions on the support of the source noise, one can recover 
the dispersion relation (the ray dynamics) in some frequency interval. A mathe- 
matical statement will be that, for r > 0, Ca.b(t) is close to the Schwartz kernel 
of f2( T ) ° n where II is a suitable pseudodifferential operator (\I/DOfor short), the 
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principal symbol of which can be explicitly computed, and 0(r) is the (semi)group 
of the (damped) wave propagation. This closeness is in general only meant in the 
L 2 sense, but can be pointwisc if the damping is strong enough. It implies that we 
can recover the dispersion relation, i.e. the classical dynamics, from the knowledge 
of all two-points correlations. 

The waves are recorded at the boundary, it implies that the main part of the 
recorded noisy fields is coming from surface waves. So we end up with another 
inverse problem: how docs the dispersion relation of the elastic surface waves deter- 
mine the structure of the Earth's interior? We will discuss the effective (dispersive) 
ray dynamics of guided surface waves. The corresponding ideal inverse spectral 
problem has been solved in [7] from which we recall the main result. 

I guess that our results can be interesting because they provide a quite general 
situation in which the source is not a white noise and for which we can still recover 
the most interesting part of the Green's function (the phase). 

Let us also mention on the technical side that, rather than using mode decom- 
position, we prefer to work directly with the dynamics; in other words, we need a 
time dependent rather than a stationary approach. Mode decomposition is often 
useful, but it is of no much help for a general system with no particular symmetry. 

For clarity, we will first discuss the nonphysical case of a first order wave equa- 
tion like the Schrodinger equation, then the case of more common wave equations 
(acoustic, elastic). 

In order to make the paper readable by a large audience, we have tried to make it 
self-contained by including sections on pscudodiffercntial operators and on random 
fields. 

Contents 

• In SectionQ] we start with a quite general setting and discuss a general formula 
for the correlation (Equation ([2])). This formula will be made more precise in 
several cases in the following Sections. 

• The goal of Section [5] is to introduce the basic examples which we will discuss 
in the paper. 

• In Scction[3l we introduce a large family of anisotropic random fields and show 
the relation between their power spectra and their Wigner measures. They are 
build from white noises using the Radonification process: our random fields 
will be the images of some white noises by some suitable linear operators. 

• Section Q] is devoted to the case of a scalar field driven by a first order in time 
differential equation. This case is quite academic, but will be used later for 
standard wave equations. We will discuss first the (nonrealistic) case where 
the source is a white noise, then the semiclassical case assuming the absence 
of mode conversion. 

• In Section [5l we discuss the case of multicomponent wave equations: we will 
discuss first the (nonrealistic) case where the source is a white noise, then the 
semiclassical case using a reduction to the case of scalar fields. 

• In Section O going back to the beginning of the story, we discuss the question 
of the correlations of codas in the framework of "quantum chaos" . 

• In Section we focus on the case of seismology and discuss the remarkable 
fact that the correlation of the surface waves is enough to image the Earth's 
crust, using a waveguide model of the crust. 
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• Finally, there are four appendices: 

— Appendix A on pseudodifferential operators, 

— Appendix B on classical dynamics, dispersion relations and generating 
functions, 

— Appendix C on Egorov Theorem for long times, 

— Appendix D on time averages versus ensemble averages, i.e. ergodicity 
of the source field. 

As can been seen in the paper, several interesting mathematical problems come 
out of our discussions; let us mention 

• Mode conversions associated with an eigenvalue going into the continuous 
spectrum (Section 17.31) . 

• Extension of the large time Egorov Theorem, with Fourier integral operators 
replacing ^PDO's, in order to be able to properly discuss the case of a source 
noise localized on a submanifold. 

• The large time behaviour of the solution of a wave equation with initially 
localised Cauchy data, like in an earthquake (Section [6]). 

Acknowledgement: I would like to thank Michel Campillo and his colleagues 
from LGIT for discussions and collaborations. I also thank the referee for his great 
help in order to improve the initial manuscript. 



1 A general formula for the correlation 

In this quite mathematical Section, we present the general context and a general formula 
for the correlation which makes quite explicit the relation with the Green's function: we will 
do that on the level of operators and on the level of more concrete functions, the integral 
(or Schwartz) kernels of the operators. Our model is quite general: it includes linear wave 
equations with some attenuation and a source noise on the right hand side of the equation. 
We will consider the following general damped wave equation: 

du i - 

— + -Hu = f 1 

at e 



1.1 Operators level 

We make the following assumptions: 

1. The parameter e > can be either fixed (say e = 1) or can be a semiclassical 
parameter going to + , like h in the Schrodinger-like case, or the inverse of 
the frequency for the wave equation. 

2. The operator iH/e is the infinitesimal generator of a strongly continuous 
group (see [36], Chap. IX) 

il(t) ="exp(~itH/e)" 
acting on some Hilbert space TL. 

3. Q(t) is a contraction group: there exists T at t > 0, called the attenuation time, 
and C > so that the operator norm of fl(t) satisfies 

Vi > 0, ||fi(t)|| < Ce-* /Tatt 
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4. The source noise f is a stationary (see Section random field f : R — > 7i. 
We assume that E(|jf(0)|| 2 ) < oojp) and the covariance of f is given by(0) 

E( |f(s) ><f(s')| ) = K(s-s') 

with K{t) G £(^0&- a continuous function of t with compact support {\t\ < 
to} (0). 

Definition 1.1 The causal solution of Equation (QP is the solution given by u = 
Gf, with G a linear operator, satisfying u{., t) = if /(., s) vanishes for s < t. 

For physical reasons, it is appropriate to take for u the causal solution of Equa- 
tion {!]). The field u is then stationary, and given by 



poo 

u(t) = / fi(s)f (t - s)ds . 
Jo 



Using Assumptions 3. and 4, one can show that u(t) is almost surely defined, and 
its covariance, defined by 

C(r)=E(|u(0)xu(-r)|) , 
is a continuous linear operator on TL. 
Theorem 1.1 The covariance C(r) is given by 

• forr > t , 

roc nt 

C(r) = n(T)n with n= / n(s)ai*( S )ds and C= / n(-u)K(u)du . 

JO J -to 

(2) 

• for t < 0, 

C(r) - (C(-r))* . 

Assuming f to be ergodic (see Section l3"3j) . the covariance C(r) is also given almost 
surely by the time average: 



C(t)= lim i [ T \u(t)><u(t-r)\dt 



In order to perform effective passive imaging, people use the previous formula for a 
large enough T and ... computers. 



1.2 Kernels level 

In order to have concrete functions and not only distributions, we will make some 
smoothness assumptions: 

• X is a smooth compact manifold of dimension d with a smooth measure \dx\. 
For applications to seismology X will be the Earth or (in Section[7]) the Earth's 
surface. 

X E denotes the ensemble average or expectation value 

2 We use Dirac's bra-ket notation, \f > is a vector, < f\ is the Hermitian scalar product with 
/, so that |/ >< g\ is the rank one linear operator: x — >< g\x > / 

3 C(E,F) is the space of continuous linear maps from the Banach space E into the Banach 
space F with the norm topology and C{E) := C(E, E) 

4 The assumption on compact support could be replaced by some suitable decay in t. It means 
that the source noise is uncorrelated for large enough times 
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• u(x, t), x G X, t € M. is the field (scalar or vector valued) with values in C N 
(or R N ). 

• The generator of the group O(t) is an elliptic (see Definition 110. ip differential 
(or pseudodifferential) linear operator H acting with a dense domain of a 
Hilbert space H C V' (X, C N ). 

• The source f is a smooth stationary random field onlxl with values in 
C N (or M. ). Its matrix valued correlation kernel is the smooth function on 
XxXxl, with values in End(C Ar ), given by 

E(f(z, s) <8> f*(y, *')) = K(s - s', x, y) . (3) 

Definition 1.2 The propagator or Green's function P is the Schwartz kernel of 
Q(t) defined by: 

(fi(t)v)(aO = f P(t,x,y)v(y)\dy\ . 
J x 

The propagator P satisfies 

/ P(t,x,y)P(s,y,z)\dy\ = P(t + s,x,z) 
Jx 

which comes from the group relation Q(i + s) = Q(t) o f2(s). The causal solution of 
Equation (fT]) is then given by 

u(x,t)= ds P(s,x,y)f(t - s,y)\dy\ . 
Jo Jx 

Due to the ellipticity of H, this is a smooth function. 
The kernel of C(r) is defined by: 

C AjB (r) = E(u(A, 0) <E> u(B, -r)*) , (4) 

where A,B £ X and CU,b(t) G End(C JV ). 

From Equation @, we get an explicit integral formula for Ca,b(t) in terms of 
P(t,x,y) and K(t,x,y); Equation Q can be written for r > i as follows^: 

C AjS (r) = [0(r)n](A, J B) (5) 

2 Examples 

/n t/iis Section, we will present the examples of "damped wave equations" which will be used 
in this paper. We start with concrete examples and go then to more general semiclassical 
examples in the scalar case as well as in the multi- component case. 

2.1 Damped Schrodinger equation 

Let X be a smooth compact boundary-less Riemannian manifold whose Laplacc- 
Beltrami operator is denoted by A. Let us give »:X^Ia smooth non negative 
function, V a smooth real valued function on X, ha non negative constant, and 
K = -hi 2 A + V(x), and take: 

h 

— {ut + a(x)u) + Ku = g . 
i 

5 If R is an operator we will denote by y) its Schwartz kernel 
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It is a particular case of Equation (TT|) where e — h, H = K + ia(x), f = &g and 
the Hilbcrt space is L 2 (X). 

Let us note for future use that, if h — > 0, the principal symbol (see Appendix A, 
Section ISTTjl of our equation reads oj+ ||£|| 2 + V(x) where our (extended) phase space 
is T*(X x R) with canonical coordinates (x,£,t,uj) and a(x) is a subprincipal term 
entering in the transport equation, but not in the classical dynamics (sec Appendix 
A). The attenuation time can be taken as T a tt = 1/ inf^gx a(x). Let us note that 
this equation has no physical meaning, it is just a kind of toy model. 

2.2 Damped wave equations 

As in the previous Section, X is a compact Ricmannian manifold. Let us start with 

u tt + au t - Au = /, (6) 

(a is a smooth function which satisfies inf xS x a{x) > 0). It corresponds to Equa- 
tion IU if we set 



and 




The Hilbert space is defined by the energy norm ||(ui, 1/2) || 2 = J x (\ui\ 2 + \u2\ 2 ) \dx\. 

Remark 2.1 It follows from ] 19$ that, if a{x) > everywhere and there exists some 
T > so that the support of a(x) crosses every geodesic arc of length > T , then 
Q(t) is a contraction group. 

One can replace the very simple wave equation by more sophisticated multi- 
component wave equations like the equation of elastic waves. The general form will 
be: 

u tt + Au t -Au = 0, (7) 

where u(x,t) 6 C N , A is an elliptic system of degree 2 and A is a differential 
operator of degree representing the attenuation. 

2.3 Pseudodifferential equations 

This example will be the subject of Section [^J 

This is a generalisation of the damped Schrodinger equation (Section 12. ip ; we 
now introduce a small parameter e which, in the case of the damped Schrodinger 
equation, is Planck's constant h. The semiclassical regime will correspond to e — ► 0, 

We assume that the dynamics is generated by a pseudodifferential operator (a 
^DO, see Appendix A). Our equation looks then like: 

u t + -H e u = f 

with 

H e = Op £ (H + sHx) . 

Here Ho is real valued and elliptic, and is the Hamiltonian of the ray dynamics, 
while ^sHi < provides the attenuation. 

In the case of the damped Schrodinger equation (Section 12. ip . we have e = H, 
H Q = ||«e|| 2 + V(x) and H x = -ia{x). 
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2.4 Pseudodifferential systems 



This example will be the subject of Section^ 
The system is of the following form 

u t + -H e u = f , with H s = Op e (H + eHi) . 

e 

Here 

• u(x,t) G 

• Hq(x, £) is an Hermitian symmetric matrix. Each eigenvalue A(a;, £) of Hq(x, £) 
gives a ray dynamics associated to a polarised wave. The polarisation bundle 
is the associated eigenbundle E\(x,£ l ) 

• Hi(x,£) is a complex matrix. Let us introduce the quadratic forms q\(x,!;) 
which is the restriction of C<sHix\x) to E\(x,£). The condition q\ < ensures 
the attenuation. 

This includes 

• The high frequency limit of the wave equations ^: we have 

)■*-(! --)■ 

A similar transform can be used for more general wave equations, like the 
elastic wave equation with N = 6. 

• An effective pseudodifferential wave equation associated with stratified media. 
They arc usually 'I'DO's with a nontrivial dispersion relation (see Section [7]). 

• Phase space dependent damping included in Hi: this is usually the case for 
seismic waves. The attenuation can be created by the small-scale physics 
which is not uniform in general. 




3 Random fields and semiclassics 

Our aim in this section is to build quite general random fields with correlation distances 
given by a small parameter e . It seems to be natural for that purpose to use e— pseudodifferential 
operators. We will see how to compute the generalised power spectrum using Wigner mea- 
sures. We will not restrict ourselves to time- dependent random fields. The subsection \3.5\ 
is devoted to the time-dependent case and we introduce there the notions of stationarity 
and ergodicity. 

3.1 Covariance of Random fields 

A random field is simply a random variable with values in some functional space, 
like a Hilbert space H. We will consider random fields f = f(uj) € TL where Q 
belongs to fl, the sample probability space. We will assume that / € L 2 (fl,H) and, 
without loss of generality, that our fields have mean 0: E(/) = 0. 

The covariance of / is the operator C = E( |/ >< f \ ) on TL. This operator is 
Hermitian and positive. TL will often be a space of functions, and we will consider 
the Schwartz kernel [C](x,y) of C. 

If TL is finite dimensional, it admits random fields for which the covariance is 
the identity operator. We call them white noises. 
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3.2 Cylindrical random variables 

For this section, see f27[ Part II] and JTEl Chap. IV] 

Definition 3.1 A cylindrical random variable / on the Hilbert space TL is given by 
the following data: for each continuous linear map A : TL — > E with AmiE < oo, 
let f a be a random variable with values in E. These random variables must satisfy 
the following compatibility condition: if I : E — > F is linear, /; .4 = I ° f A (at the 
informal level, we can write f 'a = Ao f ). 

In particular, if g : TL — > C is of the form g = F o A with A : TL — *■ E linear, 
dhaE < oo and F e L°°{E), we can define E(g) := E(F o f A ). 

Example 3.1 Let TL be L 2 (R/2ttZ), Q = with the product measure 

Il(i(<5(— 1) + <5(1))). Then f u {x) = Enez w ™ ei ™ X * s a cylindrical random variable 
on TL. f is not a random variable with values in TL, but it is a random Schwartz 
distribution in the Sobolev space H s (M./2ttZ) for s < —i. Indeed, for any lo £ £1 we 
have \\U\\ 2 HS = E„(l + H 2 ) s <cx). 

Definition 3.2 The Fourier transform w of a cylindrical random variable w is a 
function on TL defined as follows. For each x ^ 0, one can take IL X the orthogonal 
projector on the one- dimensional subspace Cx. To this projection is associated the 
random variable i«n 1 taking values in <Cx. The value of the Fourier transform 
at x is given by w(x) = E(exp(— i(wu x \x))). Informally, we can write w(x) = 
E(exp(— i(w\x))). 

To any random variable / with values in TL is associated the cylindrical random 
variable defined by fA = Ao f. If the cylindrical random variable comes from a 
random variable / with values in TL, this random variable is unique and we will say 
that our cylindrical random variable "is" a random variable. 

We have the important 

Theorem 3.1 (Sazonov |27| ) If the Fourier transform w of the cylindrical ran- 
dom variable w is continuous and A : TL — > K, is Hilbert-Schmidl[^| , then A{w) is a 
random variable on /C. The covariance of A(w) is the operator AA* . 

3.3 White noises in Hilbert spaces 

A white noise w on the Hilbert space TL is a cylindrical random variable on TL of 
mean and covariance given by the identity: 

VejGH, E((e\w)(w\f)) = (e\f) . 

Definition 3.3 The Gaussian white noise is the white noise on TL whose Fourier 
transform is cxp(— ||a;|| 2 ). 

Sazonov's Theorem 13 . 1 1 applies in particular to the Gaussian white noise. 

If the Hilbert space TL is not finite dimensional, the white noises are "random 
distributions" in the following sense: if TL = L 2 (X), we can take for TL a Sobolev 
space H s of negative order s < — dimX/2 and for A the injection from L 2 (X) into 
H S (X). So that w can be viewed as a random distribution in H S (X). 

6 A Hilbcrt-Schmidt operator A is an operator whose Schwartz kernel \A]{x,y) is in L 2 (X X X) 
and the Hilbert-Schmidt norm ||yl||[j-s of A is the I? norm of [A]. 
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3.4 An important example 

Take a self-adjoint operator H on Ti. with compact resolvent. Let us denote by Aj 
the (finite rank) spectral projection onto a compact interval /. The random field 
// = Ajw can also be defined in terms of the spectral decomposition (fj, Xj) of H. 
fj = X)a e/ a j(&)fj where the coefficients aj satisfies E(cij) = and E(aioj) = 5i t j. 

3.5 Stationary Random fields 

We will now consider time dependent random fields / : t — > f(t) g JC with t 6 1. 
Equivalcntly, our global Hilbert space is now TC = L 2 (R, JC). We can easily define an 
action of R on such a random field. For T g R, we define T./ by (T.f)(t) = f(t-T). 
We will say that 

Definition 3.4 / is stationary if there exists a group of probability preserving 
transformations on fl, denoted by {Xt, T g R}, so that T.(f(u>)) — f (Xt{oj)) ■ 

Let us consider the time dependent fields / : R — > JC as elements of a topological 
vector space B. Stationarity corresponds to the physical intuition that T.f is another 
random choice of the function / (see also Wiener's book [35]). 

If / is a stationary random field, the covariance kernel K(t,t') = E( \f(t) >< 
f(t')\ ) = K(t - t'), where K(s) g Hcrm(/C). 

Definition 3.5 A stationary random field f is crgodic if any measurable subset of 
f2 ; invariant by the transformations Xt, is of probability or 1. 

We have then the celebrated 

Theorem 3.2 (Birkhoff) If f is stationary and ergodic and F : Q — > R is inte- 
grable with respect to the probability measure, then, for almost every iv £ Q, 

hm i C F(X t (u>))dt = E(F) . 



3.6 Examples of stationary fields 

We will start with a Gaussian white noise w on L 2 (R, JC). Let us consider a Hilbert- 
Schmidt operator A : i 2 (R, JC) — > JC. Let us introduce the stationary random field 
/ defined by 

f(t) - A(t.w) . 
If the Schwartz kernel of A is a(s, x, y) we have, formally, 

f(t,x)= / Idyl / dsa(s + t,x,y)w(s,y) . 

JX JR 

Le us assume that Aw = J Wi A(s)w(s)ds with A(s) a Hilbert-Schmidt operator 
smoothly depending on s, with compact support C [— to/2, to/2]. From the formula 
in Theorem 13. II the correlation operator is 

K(t) = I A(t + s)A{s)ds (8) 
which is supported in [— to, to]. 
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Example 3.2 Stationary scalar noise on the real line: let us take K — R, F £ 
C^°(R, and the operator A given by w — > f R F(—s)w(s)ds. Then f is the 
convolution F*w. If we take the Fourier transform w.r.to time, we find for each 
frequency to £ R f(co) = F(lu)w(uj) ; so that E(|/(aj)| 2 ) = |.F| (w). The positive 
function lo £ R i— > |_F| 2 (aj) is usually called the power spectrum of the stationary 
noise f . 

Example 3.3 Brownian motions: ifK = R and w is the Gaussian white noise, the 
Brownian motion is defined by b(t) = J w(s)ds, which is in L 2 ([0,T]) for all finite 
T. Notice that b is NOT a stationary random field! 

3.7 Random fields from pseudodifferential operators 

The main goal of the present section is to build natural random fields which are non ho- 
mogeneous with small distances of correlation of the order of e — > 0. The noise is non 
homogeneous in X, but could also be non isotropic w.r. to directions. 

We will use as Hilbert-Schmidt operators some e-pseudodifferential operators. 

3.7.1 Random fields from pseudodifferential operators 

Let us take for e— dependent random fields on a manifold Z (which can be X x R) 
the image of a white noise on L 2 (Z) by an e-pseudodifferential operator A of smooth 
compactly supported symbol a(z,() £ S~^ s (see Appendix A, Definition lS.ip : this 
operator A is Hilbert-Schmidt. 

The correlation C(z, z') will then be given as the Schwartz kernel of AA* which 
is a *§>DO of principal symbol |a| 2 . We have 

C(z, z') = {2ire)- d k (z, + 0(e-^), 

where k(z, .) is the Fourier transform with respect to C of |a| 2 (z, £). 

This construction gives smooth random fields which can be localised in some very 
small domains of the manifold Z, which are non isotropic (the covariance kernel 
C{z, z') for z' close to z, depends on the direction and not only on the distance 
\z — z'\) and which have small distance of correlations of the order of e. Moreover it 
will allow to use techniques of micro-local analysis with the small parameter given 
by e. 

3.7.2 Power spectrum and Wigner measures 

The Wigner measures represent the phase space energy densities of a family of 
functions depending on a small parameter e. 

Definition 3.6 If f = (f s ) is a family of functions on Z, bounded in Lf oc (Z), the 
Wigner measures Wj of f are the measures on the phase space T* Z defined, for 
if £ C™{T*Z), by (see Appendix A, Definition ]^ 

J tpdW} := (Op e (p)/ e |/ e ) . 

The measures dWf are the phase space densities of energy of the functions (/ e ). 
We now define: 

7 If X is a smooth manifold (X, M) denotes the space of smooth functions with compact 
support 
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Definition 3.7 The power spectrum of the random field f = (/ e ) is the phase space 
measure Vj defined by: 

V) = E(W E f ) : 

the power spectrum of a random field is the average of its Wigner measures. 
Proposition 3.1 The power spectrum V of f e = Op(a)w, satisfies: 
V) = (2ire)~ d \a\ 2 (x,Z)\dxdZ\ + 0(e-( d -^) . 

Proof. - 

Let us put A = Op e (a), we have 

(Op e {(p)Aw\Aw) = (A* Op £ ((p)Aw\w) 
and VB, E((Bw\w)) = trace(B). We get 

ipdWfj = tr&ce(A* Op £ ((p)A) 

which can be evaluated using the \l/DO calculus as 

E (j tpdWf^J = {2ne)' d J tp\a\ 2 dxdi + Ofe"^) . 

□ 

3.7.3 Time dependent random fields 

If Z = X x R is the space-time, we will denote by (x,£) (resp. (t,ui)) the canonical 
local coordinates in T*X (resp.T*R). We will take our random fields as before: 
/ = L £ w; the symbol I of L e = Op E (l) is assumed to be l(x,£,w) G S~^(X x R) 
(independent of t) . The field / is stationary. 
In this case, the correlation is given by: 

K(x,y;t) = [LL*](x,y;0,t) (9) 

which is the Schwartz kernel of a ^DO of principal symbol \1\ 2 (x,£;uj). 

Lemma 3.1 The covariance operator K(t) is a ^DO in the space ^~^ s of prin- 
cipal symbol 

-!- /|Z| 2 (z,£;u,)e to / e da, . (10) 

4 Correlation for 1-component wave equation 

In this section, we will study the case of a scalar field given as in Section \2.3l First we 
describe, as a warming up, a case where the correlation is exactly computable, then we 
compute the asymptotic of the correlation in the "semiclassical case". 

4.1 An example of an exact formula for the correlation 

Let us consider the 1-component "wave equation" (with e = 1) 

du * 

— — h ku + iHqU = f 

dt 
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where k > is a constant and Ho is self-adjoint on K. with a compact resolvent. 
In particular f2(t) = e~ kt U(t) where U(t) is the unitary group U(t) = cxp(-itHo). 
Let us assume that f = $ * Aiw, i.e. f(t) = Aj (J <f>(t — s)w(s)ds) , with ?« a 
Gaussian white noise on L 2 (R,JC), $ € C£°(] - h/2, t Q /2[; E) and A/ = x/(iJ ) is 
the finite rank spectral projection of Ho associated to the compact interval /. The 
operator A : w — > f„G>(—s)Aiw(s)ds is clearly Hilbcrt-Schmidt. From Equation 
([8]) , the correlation is given by K(t) = ^{t)Ai with*(t) = / R $(t + s)$(s)ds. Using 
Theorem ll.il we get 

Theorem 4.1 We have the following simple formula for t > t a : 

C{t) = if(io)fl(r) , 

with F{lo) = xii^Ft—tuii^ (t)e kt ) pj, which is an exact relation between the corre- 
lation and the propagator P Q of the wave equation without attenuation. 

Remark 4.1 At least formally, we can take $ = 5(0), I = BL then f = w and 
F = 1. TTie correlation is then equal to the Green's function up to the factor l/2k. 

Remark 4.2 Assuming that H is "real valued" (i.e. Hf = Hf),we have Ca,b{—t) = 
C a ,b(t). 

4.2 Semiclassical assumptions 

Let us start with the semiclassical Schrodinger like equation 

i ~ 

Ut + -Hu = f 

e 

where 

• H is admissible (Dcfinition llO.ll in Appendix C): H is an e— pseudodiffcrential 
operator: 

H := Op £ (tf + £ffi) 

with 

— The principal symbol Hq(x,£) : T*X — > R, which gives the classical 
("rays") dynamics, belongs to £ m , with m > 0, and is elliptic of degree 
m (see Definition 110.11) . 

— The sub-principal symbol Hi(x,£,) belongs to S m _i and admits some 
positivity property which controls the attenuation: there exists k > 0, 
such that 

hi(x,£) := 9fHi(x,0 < -fc, V(x,£) g T*X . 

• The random field / is given by / = 0~p e (l(x, £, u>))w with w a white noise 
on X x K and with I smooth, compactly supported w.r. to (x, £) and whose 
Fourier transform w.r. to lu is compactly supported in ] — To, Tq[. The power 
spectrum of / is (27r£)~( d+1 ) |Z| 2 (a;, u>). The covariance operator Kit) is a 
^DO given by Lemma 13. II and supported in the interval ] — Tqe,Tqe[. 

• We introduce some compact energy interval I so that Support (I) C 7J ~ 1 (J). 

Lemma 4.1 Under the previous assumption on Hi, we can take for T a tt (defined 
in Section ] any number > 1/k provided that e is small enough. 

8 J 7 is the Fourier transform T(g)(u) = J* exp(— ius)g(s)ds. 
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Proof. - 



j t (v(t)\v(t))=2R(v(t)\-iH 1 (t)v(t)) 

and we use Garding inequality (see [10]): if a > 0, Op e (a) > — C for any 
C > and e small enough. 



□ 



4.3 Ehrenfest time 

Definition 4.1 (Ehrenfest time) Let X be a smooth compact manifold of dimen- 
sion d. Let us consider a smooth elliptic Hamiltonian Hq on T*X. We fix also some 
compact energy interval L . 

• The Liapounov exponent A/ = A e]0, +oo[ is any real number^) for which 
the differential of the Hamiltonian flow $ t of Hq satisfies the following uniform 
estimate: 

3C > 0, Vz e T*X with ffoOO 6 I, Vi > 0, j|d$ f (z)j| < Ce At . 

• Let its denote by e > i/ie semiclassical parameter. The Ehrenfest time, 
Tfehrcnfcst e]0, +oo] is defined as 

^Ehrenfest := | logs|/A . 

The intuition behind the definition of the Ehrenfest time is the following one: if C 
is a phase space cell of diameter ~ e, the diameter of the cell 4>t{C) remains <C 1 
for any time t smaller than Tshrenfest- We will need a smaller time: 

Definition 4.2 Let us give < 7 < h. The time T 7 = T 7 (/) is defined as 

Tf '■= ^ _ Tj ^Ehrenfest ■ 

The previous definition will be used via the 

Lemma 4.2 Assume that we have a flow $t on a compact manifold X with a 
Liapounov exponent A, then we have the following uniform estimate 

Va G N d , VA' > A, \D a $ t (z)\ = 0(e A '^) . 

4.4 Main result (N = 1) 

We get the main result: 

Theorem 4.2 With the assumptions and notations of Section \4-2\ and 7 < ^, the 
correlation is given, for r > 0, by C(t) = O(r) o Op e (7r) + R where 



n(x,0= J exph J ^(F x )($ s (x,0)ds) If (*-*(*, 0, 0) * 

belongs to S^ 00 , and is supported in Hq 1 ^). 

2 7 



'For technical reasons, not taking the infimum of these A's is more convenient 
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• 1. [No assumption on the attenuation k] we take 

k 



l = li = 



2(A + fc) ' 

and, for any a < we have the Hilbert-Schmidt estimate 

||fl|| H -s = 0{e a - d / 2 ) 

(the operator fi(r) o Op e (7r) is Hilbert-Schmidt with an Hilbert-Schmidt 
norm of the order of E~ d l 2 if ir is not = 0); 

2. [Large attenuation] if k > 2d A, we take 

fc-2dA 

7 = 72 



2(A + fc) ' 

and, for any [3 < 2j2, for (A, B, r) generic (see Definition \9.1\) , we have 
the point-wise uniform estimate 

[R}(A,B) =0(e l3 ~i) . 

Proof. - 

From Lemma 18. 11 wc know that the operator C introduced in Theorem 
11.11 belongs to ^~^ s and his principal symbol is |/| 2 (x, £, — H (x, £)). 
The correlation operator C(r) is given (Theorem II. 1|) by C(r) = f2(r) o 
J °° Q,(u)CCl*(u)du. We split the integral J °° = J" 7 + J*^°. 

• The first integral is estimated using the long time Egorov Theorem. 

• The second one is estimated using the exponential decay of fl(t) 

The main term comes from Theorem 110.21 with the formula for do 
given in Theorem 110.11 We split the remainder R into R = R± + R2 
with 



Ri = O(r) / n(s)ds 
Jo 

and ri(s) e £ 27 *l!° (using Theorem EE2) 



p 30 

r 2 = J n( s + T)cn*(s)ds 



with C G KlZ ■ 
1. The Hilbert-Schmidt estimate: from Lemma 18.21 wc get 

||i?i||H-s = 0(£- d/ V 7 |log£|) . 

Using the fact that the Hilbert-Schmidt norm of C is of order e" 
we have 

/ />oo 1 

||i? 2 ||H-s = e- d /' I e- 2ks ds 



Using the explicit value 7 = 71, we get the result. 
2. The point-wise estimate: 
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Lemma 4.3 The Sobolev H s norm of a distribution u in X is, 
uniformly in e, equivalent to 

E-||#'/ m u|| ia + ||u|| La . 

Taking s > d, using Sobolev's embedding Theorem on X x X 
and the fact that H commutes with f2(s), we get the point-wise 
estimates 

[R 2 ](A,B) = O f e-- 3 " y °° \\n( s + T)H s / m CH s/m n*(s)ds\\ li . s ds 
and 

[R 2 ](A,B) =o(e- 2s e^-^ k ' K ) . 
The result follows from the optimisation 7 = 72 ■ 

□ 

4.5 Applications: recovering the Hamiltonian from the Cor- 
relations 

Theorem 4.3 Let U be the set of (£?,£b) G T*X so that there exists t > with 
£b), — Ha(B, £b)) 7^ 0. For any tq > ; the restriction of the Hamiltonian 
Hq to U can be recovered from the knowledge of C(r) for < t < tq. 

Remark 4.3 The physical meaning of the previous assumption is that the ray dy- 
namics and its Hamiltonian generator Hq can be recovered in the domain where the 
rays cross, in the past time, the support of the power spectrum of the source noise. 

Proof. - 

Let us assume, for simplicity, that d 2 H /d£, 2 is non degenerate. For 
r > small enough, there is no conjugate points on the unique trajectory 
from B to A in time r contained in From Theorem 14.21 and 

Appendix B, we get the following expression: 

C a ,b(t) = (27Tie)~ d/2 c(A,B,T)e tS ^ B ^/ £ + R{A,B,t) 

with 

• c(A, B, t) > if (B, £b) <eU (this follows from Section EH] and the 
fact that the principal symbol w of II is > 0) 

• The L 2 (X x X) norm of R(., ., r) is 0{e a ~ d ' 2 ) with a > 0. 

The Theorem follows: S is a generating function for <3? T , the flow $ T 
determines Hq up to a constant. The constant can be extracted from 
S which is given as the action integral J £dx — Hgdt where 7 is the 
trajectory from B to A in time r. 

Without the simplifying assumption, the Theorem still holds: the 
H-S norm of the "Fourier integral operator" fl(r) o Op e (7r) is of order 

□ 
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5 The case of mult i- component wave equations 
5.1 Wave equations: source white noise 

Let us take the case of a scalar wave equation with constant damping; the closest 
results were derived in [331 [S] • We will make a formal calculation in the spirit of 
Remark |4. II We will consider the wave equation Uu + 2au t — Ait = /, with 

• a > is a constant damping coefficient 

• A a Riemannian laplacian in some Ricmannian manifold X, possibly with 
boundary: 

A = gV(x)d ij +b i (x)d i 

which is sclf-adjoint with respect to \dx\ and appropriate boundary conditions; 
in fact we could replace the Laplacian by any self-adjoint operator on X! 

• / = f(x, t) the source of the noise which will be assumed to be a scalar white 
noise (homogeneous diffuse field): 

E(f(x, S )f(y,s')) = 5( S -s')6(x-y) 



In order to get a readable expression, it is convenient to introduce Q = y/— (A + a 2 ) 
We get then easily: 



u(x, t) 



ds I e 
Jx 



sin sQ 



Q 



(x,y)f(y,t- s)\dy\ 



where sinsQ is defined from the spectral decomposition of Q (any choice of the 
square root gives the same result for sln ^Q ). We define 



G a (t,x,y) = Y(t) 



_ at sin tQ 



{x,y) 



with Y the Hcavisidc function. We will call G a the (causal) Green's function. Wc 
can rewrite: 



u{x, t) 



ds I G a (t- s,x,y)f(y,s)\dy\ 
x 



Let us assume now that / is an homogeneous white noise and compute formally 
the correlation Ca.b(t). 
We get quite easily, using 

sinasin/3 = i(cos(a — (3) — cos(a + /?)) : 



Ca,b{t) 



4a 



(Q z +a z y L cosrQ 



a sin \t\Q 



Q 



Taking the t derivative, we get the simpler formula 
d 



dT 



Ca,b(t) 



^G a {r,A,B) forr>0 
±G a (-T,A,B) forr<0 



(A,B) 



(11) 



(12) 



This "exact" relation between the correlation and the Green's function in the case 
of a white noise source field is the starting point of many works in seismology. It 
has been derived in various ways by several authors in particular R. Weaver and O. 

Lobkis in EU [M nn. 
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Remark 5.1 The previous calculation can be extended to the more general wave 
equation 

u tt + ku t ~ Am = / 

with k > constant, A self-adjoint and f a white noise. In particular, Ca,b(t) is 
an even function of r . This is the basis of 

• Clock synchronisation \22§ : if the clocks at the stations A and B are not 
synchronised, the correlation is even with respect to a shifted origin of times. 

• Test of applicability of the theory: in the semiclassical regime studied later, 
approximate evenness holds. In fact the oscillating part is even because the 
corrections to the Green's function do not affects the phases. In Theorem \4-.S\ 
the principal symbol 7r is > 0. This is basically due to the time reversibility of 
the ray dynamics. 

5.2 General wave equations: semiclassics 
5.2.1 Assumptions 

We want to derive results similar to those of Theorem 14.21 in the case of (multi- 
component) wave equations. For that, we will introduce a matrix version of what is 
done in Section [4] and use the results of Section [8?7l Let us start with the following 
case of Equation |T]) : 

-u t +Hu=-f (13) 

I i 

where 

• H is a matrix of e— pseudodiffcrcntial operators: 

#:=Op £ (ff +£#i) 
with (see Appendix A, Section UTT)) 

— The principal symbol H (x, £) : T*X -> Herm(C Ar ) belongs to X^f^™ 
and is elliptic of degree m with the eigenvalues Ai (x, £) < • • ■ < Xj (x, £) < 

■ • ■ < X n {x, £) of constant multiplicities mj; the corresponding cigenspaces 
Ej(x,£ l ) C <C N , of dimension rrij, are called the polarisations. This as- 
sumption is satisfied for the elastic wave equation in a manifold of dimen- 
sion 3 with N = 6: 2 eigenvalues ±As of multiplicity 2 (the S— waves) 
and 2 eigenvalues ±Ap of multiplicity 1 (the P— waves). 

— The sub-principal symbol H\{x,C) £ EjjJ, ' C admits some positivity 
property which controls the attenuation: there exists k > 0, such that 
the sub-principal symbols H\.j, j = 1, • • • , n, satisfy 

< -k < . 

• The random field f is given by f = Op e (l(x, £, w))w: 

— The noises w = (w\, ■ ■ ■ , Wjsr) arc independent white noises onlxl 

— The matrix I is smooth, compactly supported w.r. to (x, £) and its 
Fourier transform w.r. to u) is compactly supported. The power spectrum 
of / is (2Tre)~( d+1 hl*(x,Z,u). 
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5.2.2 The main result 



Using the gauge transform reducing H to the block diagonal normal form of Section 
18.7.21 we can decompose the correlation into blocks C\ b (t), j, k = 1, • • • , n. 

Theorem 5.1 With the assumptions of Section \5.2.1\ the correlations C^ B (r) are 
given 

• if j ' = k) as * ra Theorem \4-^ with Hq = Xj, 

• if j ^ k and l(x, £,w) vanishes for x near B: C^ B (r) = 0(s°°). 

The proof makes a strong use of the normal form of Section 18.7.21 we use a block 
decomposition of C = (C^) in Equation (J5J) and the block diagonal reduction fi(£) = 
(Oj(i)). We are reduced to the study of 

/>oo 

Jo 

If i = j, we are in a situation very close to that of a scalar Hamiltonian (Section 
14. 4|) . while if i ^ j, we can use Lemma \l 0.2 1 

In particular, for the elastic wave equation, if the source noise is supported away 

• PS 

from B, the correlation C A ' B between S-waves and P-waves is very small. 



6 Coda's correlations 

In [5], M. Campillo and A. Paul computed the correlation of coda waves and saw 
the emergence of the Green's function. Let us try, rather informally, to discuss this 
observation, even if we do not have mathematical results on it. 

Let u{x,t) be the solution of the wave equation with initially localised Cauchy 
data (delta functions, no randomness in the Cauchy data) at time t = (an " earth- 
quake") The coda is the field u(x,t) in some window [Ti, T 2 ] with T\, T 2 >> 0. Let 
us look at the correlation 

Ca,b(t) = lim — -i— f 2 u(A,t + r)u(B,t)dt 
ri,r 2 -»oo J 2 — J-i Jti 

assumed to exist. Introducing the propagator P, we get: 

fT-2 \ 



CaAt) = I \dx\P(r,A,x) [ lim — ^— f 



u(x, t)u(B, t)dt 
J 



So that the question is reduced to the study of the limits of T2 ^_ Tl u(A, t)u(B, t)dt 
as Ti,T-2,T2 — T\ — > oo. If we have some nice (smooth?) limit K(A,B), then 
Ca,b(t) = [n(T)oK](A,B). 

Looking at some semiclassical limit, this is closely related to the following result 
by R. Schubert |26j : if the classical dynamics is hyperbolic and the time t of the 
order of the Ehrenfest time, the smoothed Wigner measures ("Husimi" measures) of 
u(t, .) are equidistributed on each energy shell. Does this uniformity of the energy 
distribution extend to the wave functions themselves, in particular to the phases? 



7 Using surface waves in order to image the inner 
crust 

In real applications to seismology, the fields are recorded at the surface of the Earth. More- 
over the main contribution to the correlations comes from the surface waves: 
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• They have a smaller time decay than the body waves. 

• The source of the noise itself (interaction with atmosphere and ocean) is located at 
the surface. 

The Earth's crust acts as a wave guide. This implies that we really have to look at a 
2D version of the wave equation: the wave equation for surface wave is an effective wave 
equation on the surface. This equation is given in terms of eigenvalues of a vertical Sturm- 
Liouville equation. Solving an inverse spectral problem for these Sturm- Liouville equation 
allows (in principle!) to image the crust. The inverse spectral problem is solved in [7]?. 

In this Section, we will discuss these effective Hamiltonians in the scalar model of the 
acoustic wave equation and show how one can reduce the problem of imaging the Earth's 
crust to an inverse spectral problem. 



7.1 A mathematical model 

We work locally in X = {(x, z) £ R 2 x E | z < 0} so that z = is the surface. We 
will consider the very simple case of an acoustic wave equation near the origin of 
X . Let us give a function 

AT(x, Z) : R 2 x K_ -> R + 
which satisfies the following assumptions: 

1. There exists > and Z < 0, so that N(x, Z) = for Z < Z . 

2. For every x 6 R 2 , M z N(x, Z) = N (x) > and N (x) < for x close to 


3. N is C 1 w.r. to Z and the map x — ► JV(x, .), from K 2 into the Banach space 
^([ZcO^R), is smooth. 

Let us introduce a small parameter e > and the following acoustic wave equa- 
tion where 

n(x,z) = 7V(x,-) : (14) 



e 



Utt — div(n gradu) = 

§!(x,o) = o 



(15) 



We plan to see that Equation (fT5|) admits, as e — > 0, asymptotic solutions 
of frequency of order e _1 located near the surface. Moreover these solutions are 
determined by solving an effective pseudodifferential equation on the surface dX = 
R 2 x {0}. In order to show that we can apply Theorem 18.31 of Appendix A, we will 
first introduce the operator Hq in Section [7. 31 It will be a Sturm-Liouville operator 
in the variable Z. 



7.2 Discussion of the assumptions on the model 

The scalar acoustic wave equation (|15j) is a simplified model of the elastic wave equa- 
tion. Wc could work with the clastic wave equation, but it will be more complicated 
and the main ideas would be the same. 

The function n represent the square of the local propagation speed (group veloc- 
ity). The formula (fT4)l means that the speed admits fast variations in the vertical 
directions and is smooth in the horizontal direction. It is a plausible assumption 
that the the geological structure near the surface of the earth is "stratified" . 

Assumption (1) on A^ means that the speed is constant far from the surface. It 
is an assumption which could be removed at the price of a technical work. 
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Assumption (2) implies that the speed is somewhere smaller than the speed at 
large depth. This kind of assumption is usually satisfied in seismology where the 
speed of elastic waves into sediments layers is smaller than the speed inside the 
rocks below the sediments. 

Assumption (3) means that we need less regularity in the vertical direction. It 
could be relaxed using piccewise smoothness w.r. to Z, but at the price of a rather 
technical work: the domain of the Sturm-Liouville operator introduced below would 
then depend on x. 

7.3 A family of Sturm-Liouville operators 

Let us consider, for each (x, £) <G T*dX, the self-adjoint differential operator L x ,t 
on the half line Z < 0, with Neumann boundary condition at Z = 0, defined by: 

L^v := - A (V(x, + iV(x, Z)\tfv . (16) 

The domain D C L 2 (R-) of L x>c is 

D := H 2 (E.~)D{v | v'(0) = 0] 

which is independent of (x, £). 

The spectrum of L X; £ consists of a finite discrete spectrum and a continuous 
spectrum [iVoo|£| 2 , +oo[. i x ,4 admits, for £ large enough, a non empty discrete 
spectrum of simple eigenvalues 

A^ (x)|e| 2 < Ai(x,0 < • •• < < ■ ■ ■ < A fc (x,0 < A^l 2 , 

which depend smoothly of (x, £). In order to see that, we can interpret Cn = 
|£|~ 2 L X .£ as a semiclassical Schrodinger type operator with an effective Planck con- 
stant h = |^| _1 and a principal symbol 

Px (Z,C)=N(k, Z)(C 2 + 1) 

which admits a "well" because of Assumption 2 in Section [77T1 on N. 

We should however take care of the fact that the number k of eigenvalues depends 
on (x, £) and goes to oo as £ does. This will give birth to mode conversions between 
surface waves and body waves. To my knowledge, such mode conversions have not 
yet been studied from a mathematical point of view. 

7.4 Wave guides as an adiabatic limit 

We will rewrite Equation (fT5|) as an operator valued pseudodiffcrcntial equation as 
in Appendix A, Section [7^1 in the spirit of the book [55], Chapter 3 on "space- 
adiabatic perturbation theory" . 
Putting 

u(t, x, z) — v(t, x, -) 

e 

we get the following wave equation for v: 

e 2 v tt - (d z (Nd z v) + e 2 NA x v + e 2 grad x A^.grad x ij) = , 

which can be rewritten, using Weyl quantisation in the variable x: 

e 2 v u + Op e (Z x ,£ - ie£.grad x Af ) v = . 

Putting ffo((x, £)) = Lx,£ anc ^ ^i(( x jC) = — *C-g ra d x Af, we will apply Theorem 
18.31 in Appendix A. Let us choose, for (x, £) close to 0, some eigenvalue Xj(x, £) of 
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L x ^ with a normalised cigenfunction v?j(x, £, Z). Then wc have the existence of a 
propagation mode of the form 

v(t,x,z) = w(t,x) (tpj(x,£,z/e) + £•••) 

where w satisfying the scalar wave equation: 

e 2 w tt + Op £ (\j(x,£) + ■ ■ -)w = . 

For each eigenvalue Xj , we can then apply the results of Section [5] 

From the passive imaging method, we can hope to recover the effective Hamilto- 

nians ±^/Aj. Solving an inverse spectral problem for the Sturm-Liouville operator 
allows in principle to recover the function n(x, z) which encodes the physics of 

the propagation medium. 

From the paper [7], we know that the function N can be recovered from the 

large £ asymptotic of the A/s. 

8 Appendix A: a review about pseudodifferential 
operators 

We review the basic definitions and properties of pseudodifferential operators (ty DO's). 
Pseudodifferential operators (first called singular integral operators by Calderon and Zyg- 
mund) where developed in order to get a class of operators generalising linear differential 
operators and containing approximate inverses of elliptic operators. They became a ba- 
sic tool in the theoretical study of linear differential operators with variable coefficients. 
A more general theory was then developed including a small parameter, here denoted by 
e €.]0,£q\. Several textbooks are now available Wl [771 \30f . 

8.1 Basic calculus 

We will define the pseudodifferential operators ('FDO's) on R d . fDO's on manifolds 
arc defined locally by the same formulae. 

Definition 8.1 (Classical symbols) • The space of symbols of degree k 
is the space of smooth functions p : T*M. d — > C which satisfy 

Va,/?GN d , \D2Dl P (x,t)\<C atfl {$ k , 

with (0 = i + U\\. 

• A classical symbol of degree m is a family of functions 

p e : T*R d -> C 

depending of some small parameter e > which admits, as e — > 0, an asymp- 
totic expansion 

oc 

Pe = ^Vpj(a;,£) 

with pj <E S m . We will denote this space by 5™ ass . 

Definition 8.2 (Pseudodifferential operators) A classical e-pseudodifferential 
operator P £ (a ^DO) of degree m on M. d is given by the (distribution) kernel 

[p e ]( X ,x') = pke)-* J e«*-* 'w Pe ( £ ±^,e) m 
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where p £ (x, £), the so-called (total) symbol of P, is in S™ ass . 

We will denote P e = Op E (p e ) and by ^"lass the space of classical ^DO's of degree 
m. The operator P £ is called the Wcyl-quantisation of p £ . 

The operators P e act from into C°° . The kernel of P e is given by: 

[P £ ](x,x>) = (2ne)- d p e (i±^, (17) 

with p £ the partial Fourier transform of p £ (x,^) w.r. to £. Very often, one is only 
able to compute the symbol po which is called the principal symbol of P. The term 
Pi is called the sub-principal symbol. 

The most basic fact about $DO's is the fact they form a (non-commutative) 
algebra: if P e = Op £ (p £ ) and Q £ = Op £ (q £ ), we have P £ Q £ = Op £ {p e q e + 0(e)). 
The composition formula for the total symbol is called the Moyal ^-product Op e (a)o 
Op (b) = Op e (a * b) and is given by 

oo 

a*b = ab+ ^{a, b} + e 3 P-j(a, b) 

where 

• {a, b} is the Poisson bracket 

• Pj(a,b) is a bi-linear bi-differential operator, homogeneous of degree j with 
respect to a and b. 



8.2 The symbols 

In order to formulate the Egorov Theorem for long times, we will need more sophis- 
ticated classes of symbols. 

Definition 8.3 (S™ ) For any real number m and any 5 with < S < ^, a symbol 
a £ S™ is a smooth function on T*X depending on e which satisfies 

Vq £ N 2d , 3C a > 0, V(x,£) e T*X 7 \D^a(x,^\ < C a e~ 5W <O m ■ 

We have S% ass C Sff*. 

We can associate to a £ S™ a pseudodifferential operator Op £ (a) £ \P™ using 
formula p7|) . Such operators obey a nice calculus, see \TU\ I14j: if, for j = 1,2, 
<2j € 5^ 3 , Aj := Op e (cij) and J := max(<5i, 82), we have A\ o A 2 = Op e (ai * a 2 ) 
with ai * a 2 € 5 , ™i+™2 an( j ac j m it s the asymptotic expansion given by the Moyal 
★-product. We have P J {a 1 ,a 2 ) £ e~^ Sl+s ^ Sp +m \ In particular [A u A 2 ] = 
f Op e ({01, Oa}) + Op e (6) with b £ £ 2(l-5 1 -<5 2 ) 5 ™ 1 +m 2> 

Theorem 8.1 (Calderon-Vaillancourt) The operators of^g extend to operators 
from L 2 to L 2 , uniformly continuous in e. The operators in 'J™ extends to operators 
from L 2 into the Sobolev space H~ m with a norm of order 0(e m ). 



8.3 Asymptotic expansions 

The basic tool in order to find a ^DO is an inductive construction based on the 
realisation of any asymptotic expansion: given a sequence aj £ S™' and li = < 
1% < • ■ ■ < lj < ■ ■ ■ and lj — > +00 as j — > 00, there exists a £ S™ such that 



E 

3=1 



e lj a-i 
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meaning that 

71-1 

Vn, a - JV'aj € el " S ? ■ 

5=1 

Such an a is uniquely defined modulo a remainder term r which lies in e°°S' s n . The 
corresponding ^DO Op £ (a) is well defined modulo 0(e°°). This implies, if m < 
and a compactly supported in x, that the L 2 — > L 2 norm of Op e (r) is 0(e°°). 

8.4 Microlocalisation 

A family (u £ ) s ^q of functions on X is said to be admissible if the L? norm on any 
compact set is at most of polynomial growth with respect to 

The family of functions (u e ) is said to be 0(e°°) near a point (x, £) € if 3</j £ 
C™(X), tp(x ) ^ and the e-Fourier transform f £ (^ £ )(^) := (27re)- d / 2 / e^l^^^dx 
satisfies .F e (v5it e )(£) = 0(£°°) uniformly for £ close to £ - 

The micro-support MS(u e ) of u e is the set of all (x, £)'s at which u £ is not 0(e°°). 
For any fDO A, we have MS(Au e ) C MS(u £ ). 

If U is an open subset of T*X and if w e and v e are admissible, we will shortly 
denote u £ = v £ + 0{e°°) in U if MS(u - v) H U = 0. 

8.5 Functional calculus and and an useful Lemma 

A simple case of the functional calculus for f DO's is as follows: 

Theorem 8.2 Let f £ C^°(R) and A = Op e (a +£aH ) £ *™ aSB with a a elliptic 

(see Section al 0.1\) and real. Then one can define f(A) using the spectral theory and 
f(A) £ ^Zj™ with principal symbol /(a ). 

Remark 8.1 The previous Theorem is usually stated for A self-adjoint (see \28\j 
Chap. 8). It can be extended to the case where the sub-principal symbol is not real, 
with a very close proof. 

Lemma 8.1 Let us consider the operator C, given by 

£:= i/ / n ( -< ) e<to/E OpZ(L(x,Z,u))dudb 

with L £ S~^! S {X x R) with L(x,t;,Lo) having an u-Fourier transform compactly 
supported as a function of (x, £, i). Here Q(t) = cxp(— itH/e) where H satisfies the 
assumptions of Section \4-2\ Then C belongs to ^~(^ s and its principal symbol is 
given by l (x,g) = L(x,£,-H (x,g)). 

The integral giving C is supported w.r. to t in the interval [—to, to]- Using a 
change of variable t — et' in the Definition of £, the proof is a simple corollary 
of the functional calculus of WO's: if A = Op e (a) with a compactly supported, 
the operator B = exp(it' H)A belongs to Vf^^ss anc ^ *he principal symbol of B is 
b = cxp(it' H())a; because the function exp(ii-) is not compactly supported, we need 
first to rewrite A = F{H)A with F £ and then use the functional calculus with 
the function / = Fexp(ii-). 

8.6 Hilbert-Schmidt estimate 

Lemma 8.2 // P £ ^g°°(X) with X a compact manifold, P is Hilbert-Schmidt 
with ||P|| H -s = 0{s- d ' 2 ). 

Proof. - 
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P*P G 00 is trace class (it has a smooth kernel) and its trace is 
easily seen to be 0(e~ d ). 



□ 



8.7 Multicomponent operators 
8.7.1 Systems of *DO 

Definition 8.4 Let E and F be 2 Banach spaces S„ F is the space of symbols of 
degree m with values in the Banach space C(E, F) of continuous linear maps from 
E into F. 

Similarly, we can define the spaces S™^ S ' F , ^™^ S ' F , ••• ■ The operators of ^™^ S ' F 
send functions on X with values in E into functions on X with values in F. 

Definition 8.5 Hq G S^'^ is elliptic if Ho(x,£) is invertible for large £'s and the 

■ ■ ^F,E 

inverse is in 2j_ m . 

Matrix valued fDO's have properties similar to scalar fDO's. 



8.7.2 Reduction 

Let H be an Hilbcrt space and D C Ji a dense sub-space. 

Definition 8.6 A matrix valued Hamiltonian Hq : T*X — > Herm(W) is an element 
o/S^ ,w : for every (a:, £) G T*X, Hq(x,£) is a self-adjoint operator with domain 
DcH. 

Let us consider a classical symbol H e = Hq + eH\ + ■ ■ ■ G S"^ S ' H with Hq a matrix 
values Hamiltonian. Let us consider near some point zq G T*X an eigenvalue \q(z) 
of constant multiplicity k and cigenspace E(z) (the polarisation space) and assume 
that Xq(z) is an isolated point in the spectrum of Hq{z). 
Then we have, inspired from [TBI I29j: 

Theorem 8.3 In some neighbourhood of zq, there exists 

• A C(H)~valued ^>DO U = Op e (£/ + eU\ H ) with U Q (z) unitary. 

• A classical symbol in S™^ s ' C X(z) = Ao(z)Id + eXi(z) + • • • 

• A symbol K with values in C(T>f)F, F) with F = E with Aq(z) ^ Spcctrum(Xo(z)), 
so that 

near zq. 

Moreover Ai(z) = ttHitt + R\ with it the orthogonal projector onto E and R\ 
self-adjoint. 

The proof is done by induction on the powers of e: at the first step, we choose an 
unitary gauge transform U(z) G C(Ti) so that 

with Aold — Kq invertible. At the next step, we try to find 

-('.!) 
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so that 

< 1 - B )*(( T /• )-*.)* ,1+£B, ^( A ° Id +£Al + ' +0(£Z) 
Collecting terms in e , we get the following equations: 

6(A Id - ifo) = 0, (A Id - Jv )c = - 7 

with 

V 7 o 

These equations can be solved with the assumption that Aold — Kq is invertiblc. 

Remark 8.2 In the Theorem \8.3[ we do not need that Op E (H) is self-adjoint (only 
0~p e (Hq) is self- adjoint), neither that the symbol of Op £ (if) takes values in bounded 
operators on Tt. 



Remark 8.3 The Theorem \ 8.3\ will be used as well to reduce a multi- component 
wave equation to scalar equation in Section \5.S\ and as tool for adiabatic limits in 
Section \7.4\ 

Corollary 8.1 Assuming that all eigenvalues Ai < • • • < Xj < • • • < A„ of Ho(x,£) 
are of constant multiplicities rrij, j = 1, • ■ ■ , n, the Hamiltonian H admits a micro- 
local normal form with diagonal blocks of sizes rrij X rrij and principal symbols of 
the j—th block Xjldc m i ■ 

The sub-principal symbol of the j—th block is H\j := TTjHiitj + Rj where Rj is 
self-adjoint while 7Tj is the orthogonal projection on Ej ; it means that, as quadratic 
forms, Hij is the restriction of Hi to Ej. 

The formula for the sub-principal symbols is a direct consequence of [13j . 

8.7.3 Wave equations 

Let us consider the wave equation ([7]) 

u tt + Au t — Au = 

with e 2 A = Op e (Ho) (Hq elliptic) and assume w.l.o.g. that —A is > so that 
A i = eV— A is also a *DO. Putting 



U := 



Aim 

2 

—ieut 



Equation {7J reads as: 



which is of the form 



e / O -Ai 
^ + ( -A, -ieA ) L/ = ' 



i 



-U t + HU 

i 



with H an elliptic ^DO. 
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8.7.4 Extension of Lemma [87T1 

Let us consider a <G S(M. k3 ,L(C )) and Hq an N x JV Hcrmitian matrix, we define 
a(H ) as 

a(J3o) := ^- / e ltH °Ta{t)dt . 
The following Lemma is proved the same way as Lemma 18. II 
Lemma 8.3 Lei us consider the operator C, given by 

£:= i/ / n (- t y tuj/£ °p x e( L ( x ^^)) dujdt 

loii/i L <E S' c i ass ' ' x R) iirci/i L(x,£,lj) having an 10 Fourier transform com- 

pactly supported as a function of (x,£,i). -ffere f2(i) = cxp(— itH/e) where H 

satisfies the assumptions of Section 1 5. ,21 Then C belongs to ^ ' and his 
principal symbol is given by Iq(x,^) = L(x, £, — H (x, £)). 

9 Appendix B: asymptotics of the Green's func- 
tion and classical dynamics 

9.1 The scalar case 

We want to describe quite explicitly the propagator P(t,x,y) = [Cl(t)](x, y) with 
Q,{£) = exp(— itH/e), H as in !4.2l We will do that in the energy interval /; it means 
that we will only describe the action of fl(t) on functions whose micro-supports are 
in Ho 1 (I). 

To the Hamiltonian Hq : T*X — > M, we associate the ray dynamics defined by 



dxj dH d£j dH 
dt d^j ' dt dxj 



(18) 



Let us denote by 4>t the flow of the Hamiltonian vector field Xh defined in 
Equation (|18| . The main result is that P(t,.,.) is a Fourier Integral Operator 
associated to the flow <p t . For simplicity, we will only describe the so-called generic 
case where P admits a WKB expansion. 

If x$,yo € X and 0t(2/O)??o) = (^OjCoX we say that Xo and yo are not conjugate 
along the ray 7(s) = 4> s (yo,r)o), < s < i, if the differential Sn — > d(f> t (0,Sr]) is 
invertible at the point (yo^o)- This implies that for (x,y) close enough to (xo,yo) 
there exists an unique (y,r)) close to (yoj%) with <pt{y,rj) = (x,£). In this case 
the action S-y(t,x,y) is defined by S 1 (t,x,y) — J^^dx — Ho(x,£)dt. S 7 (t, ., .) is a 
generating function of 4>t near (j/q, t?o), meaning that <fit(y, —dS 7 /dy) = (x, dS^/dx) 
for (y,r?) close to (y ,?7o)- 

The action S* 7 is useful in order to describe a WKB expansion of the propagator 
P(t,x,y) for (x, y) close to (xo,yo) : P(t,x,y) is a sum of contributions P 7 (i, x, y) 
of rays going from y to x in the time t. 

Assuming that x and y are not conjugate along 7, we have 

P 7 (t,x,y) = (2^)- d / 2 ff>X T (i,x,y)j . 
Moreover ao, 7 (t, x, y) > if t is small enough. 
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Definition 9.1 We will say that (x,y,t) is generic if x and y are not conjugate 
along any 7 so that 7(0) = x, j(t) = y and 7 C Hq (I). 

If (x, y, t) is generic, there is only a finite number of rays from y to x in time t 
included in Hq (I). 

Dispersive rays: we say that the dynamics generated by Hq is dispersive if (Hq)^ 
is invertible. This means that the map £ 1— > v(£) := OHq/O^ (v is the speed) is a 
local diffeomorphism. If this is the case, (t, x, y) is generic for x close enough to y. 

The dispersion relation: the principal symbol of the wave equation is w + 
Ho(x,£) : T*(X x R) — > K. It is clear enough that the ray dynamics is deter- 
mined by the hypcrsurface D := {(x, £, ui)\lu + Ho(x,£) = 0}. This surface is called 
the dispersion relation. This terminology does not mean that Hq is dispersive! 

9.2 The matrix case 

Let us assume that Hq is matrix valued as in Section 15.21 The dispersion relation 
is then the hypersurface 

det(wld- H (x,£)) = B(x,Z,w) = , 

we can (at least outside a mode conversion set), rewrite 

B(x,^)=llf =0 (cj-X j (x,0) di 
which gives M ray dynamics Xj associated to different polarisations 

Ej =her(Ho(x,£)-\ j {x,€)U) . 

10 Appendix C: long time Egorov Theorem 

The main reference for this part is [3] where the authors get only an operator norm 
estimates for the remainder. The proofs of reference [3] can be extended in order 

• to get a point-wise result under some cllipticity assumption 

• to cover the case of manifolds 

• to use a local Liapounov exponent: the Liapounov exponent in some energy 
interval / 

• to allow an attenuation term. 
10.1 Egorov theorem 

Definition 10.1 An admissible Hamiltonian of degree m > is H := Op e (i?o + 
sHi) with 

• Hq £ S m a real valued function which is elliptic 

3c, C>0, V(x, g T*X, H (x, > C(0 m - c . 

• Hi G S m _i and hi = < — k < everywhere in T*X . 

Example 10.1 A typical admissible Hamiltonian is the damped Schrodinger oper- 
ator H = — e 2 A + V(x) — isa(x) with V(x) > and a(x) > 0. 
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We define the dynamics, for t > 0, by: 

n(t) = exp(-itH/e) . 



Theorem 10.1 (Egorov Theorem) Let H = Op e (H) be an admissible Hamilto- 
nian and A(0) = Op £ (a) with a G C^ o (iJ _1 (/)) with I a compact interval. Then the 
operator A(t) = n(t)A(0)n* (t) is, for t fixed, a pseudodifferential operator whose 
symbol aOt) = ^ e J aj (t) belongs to S~^ s and satisfies 

a (t)(z) = exp ^2 J s (z))ds^j o($_ t («)) . (19) 

Let us remark that the previous result is a priory uniform only in a fixed interval 
t G [0, To] independent of e. 

Let us sketch the proof. A{t) is the value for s — t of the solution B{s) of the 
ODE 

^-Q,(t - s)B(s)Q*(t - s) = 0, B(0) = A(0) . 
ds 

We have 

■^nOt - s)B(s)n*(t - s) = n(t -s)( B'(s) + -(HB(s) - B(s)H*)) fi*(i - s) . 
ds \ e ) 

Trying B(s) = Op e (6 (s) + ebi(s) H ) and using the Moyal product, we have 

to solve 

b'(s) + -({H + eHt) * b(s) - b(s) * (Hq + eHi)) = . 
We get the following equations by looking at the term: 

— bj + {H ,bj}i - 2h\bj = c , 

with Co = and, for j > 1, Cj = X)fe=o Pk,j(pk) where the Pkj's are linear differential 
operators constructed by using the Moyal product with Ho and H\. If B(s) = 
Op e (6o(s) + ebi(s) + •••), we have 

—not - s)B(s)n*(t - s) = not - s)R(s)n*(t - s) 

ds 

with R(s) G e°°*^. 

It remains to prove that fi(t— s)R(s)n*(t — s) is 0(e°°) as an operator from the 
Sobolev space H~ s into the space H s , for any s. From the Calderon-Vaillancourt 
TheoremO it is enough to prove that L 2 -> I? norm of H s / m n(t - s)R(s)n k (t - 
s)(H*) s ' m is 0(e°°). This is clear because the powers of H (resp. of H*) commute 
with n(t-s) (resp. Q*(t-s) ). 

10.2 Long time Egorov Theorem 

We want to extend Egorov Theorem to times going to infinity when e — * 0. From 
Lemma l4~2l we get (see [3]): 

Lemma 10.1 Let a Q G C™(T*X) and I = i? (Support (a ))- Let T 7 (I) be the 

time defined in Definition \4-^\ Then, the function ao(t) belongs to S^°° with S = 

i — 7 < i for all times t G [0, T 7 (7)], and the constants in the estimates of the 

derivatives are uniform in t. Similarly, the symbolds of higher order e J aj belong to 
^■(1-25)5.-00 _ 



29 



The main result is: 



Theorem 10.2 (Long time Egorov Theorem) Let H = Op e (H) be an admis- 
sible Hamiltonian fSection \10.1\) and a £ C^°(T*X). For any t £ [0, T 7 ] and with 
uniform symbol estimates w.r. to t, the operator A{t) := fl(t)AQ*(t) is a pseu- 
dodifferential operator, the total symbol of which is given as in the classical Egorov 
Theorem \10.1\ by 

A{t) = Ve (f^ •'„,:/ >j , 

where e'aj(t) € e'^-^Sg 00 . 

More precisely, for any M and k, there exists N so that if -Rzv(i) = A(t) — 

o Pe (e; lo £? ^ e C k norm of the kernel of Rjsr(t) satisfies, uniformly for 

t G [0, T 7 ], the bound 

\\[RN(t)]\\c*{xxx) =0(e M ) . 

Remark 10.1 The long time Egorov Theorem was first proved by Bambusi, Graffi 
and Paul J7J/ and improved by Bouzouina and Robert 

Our estimate of the remainder term is better than the one given in which is 
only an L 2 operator norm. We need the ellipticity of Hq . 

Using Lemma llO.il the proof is the same as for a fixed time interval (see Section 



10.3 Integrals with two dynamics 

Lemma 10.2 Let us assume that we have two admissible Hamiltonians H± = 
Op e (-ffo.± + £-Hi,±) with k± > 0. Let us define the operator 

/>oo 

J := n+(t + T)Op £ {a)n*_(t)dt, 
Jo 

with a € C£°(T*X). If we assume that, for z £ Support(a), Ho t +(z) ^ Hq^(z) (no 
mode conversions occur in the support of a ), then 

J = sQ + (t)K + R , 

with K a ^DO of principal symbol fco = —ia/{Hv.+ — #o,-) an d R € e 00 ^ -00 . 
In particular, if a vanishes in some neighbourhood of B, [J](A, B) = 0(e°°). 



Proof. - 



There exists some *DO K so that Op e (a) = —i(H + K - KH_)* + r 
with r € £ 00 \D' -00 . J can be rewritten as 

j = e — (n+(t + T)Kn*_(t))dt+ / (n + (t + T) r n*_(t))dt . 

Integration by part gives the result. 

□ 
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11 Appendix D: stationary Gaussian fields are er- 
godic 

This section is inspired from [2J. Let us define 

1 f T 

C t {t,A,B):= — u(A,t) ® u(B ,t - r)dt , 
T Jo 

where u is the causal solution of Equation([T]). We then have 

Theorem 11.1 Let us assume that f(.,t) = Kw(.,t) where w(x,t) is a (Gaussian) 
white noise on X x M and K : 7i — > 7i. Then 

1. If K is Hilbert- Schmidt, Ct{t, A, B) — > E(C(t,A,B)), as T — > oo, almost 
surely for almost every pair (A, B) . 

2. If K is smoothing and H elliptic, Ct{t, A, B) — > E(C(r, A, B)), as T — > oo, 
almost surely in the C°°(X x X) topology. 



Proof. 



1. We will assume w.l.o.g. that u is scalar valued and real and use the 
following classical identity: if / = (fi, fi, fa, f4) : — ► M 4 is a Gaussian 
random variable, 

nhfafah) = E(/i/ 2 )E(/ 3 / 4 ) + E(/ 1 / 3 )E(/ 2 / 4 ) + E(/ 1 / 4 )E(/ 2 / 3 ) . 

(20) 

We have 

e (I|Ct(t)|||_ s ) = ^ / d< / dt ' I dAdBE(u{A,t)u(B,t-T)u(A,t')u(B,t' 

1 Jo Jo JxxX 

We apply the identity ([20|) and get: 

E(||C5r(r)|||_ s ) = ||E(C T (r))|||_ s +^ f dt f dt' [ dAdB(II+III) , 

1 " Jo JO JXxX 

with 

// = E(u(A, t)u{A, t'))E{u{B, t - t)u(B, t' - r)) 

and 

III = E(u(A, t)u{B, t! - t))E(u(A, t - t)u{B, t')) . 
For example, we have 



OO /-OO 



E(u(A,t)u(B,t'-T)) = [ ds / ds'n(s)KK*fl*(s')6(t-s = t'-s'-T)](A,B) 



whose L? norm as a function of (A, B) is estimated as 0(exp(— \t — t' + 
t|)). Using Cauchy-Schwarz inequality and the fact that a product of 
two Hilbert-Schmidt operators is trace class, we get 

|i f dt f dt' [ dAdB(II + III)\=0(l/T) . 
T Jo Jo JxxX 

Standards argument involving Tchcbichcv inequality allow to conclude. 

2. The second assertion is proved in a similar way using the Sobolcv 
embeddings, like at the end of the proof of Thcorcm llO.il 

□ 



31 



References 



D. Bambusi, S. Graffi & T. Paul. Long time semiclassical approximation of 
quantum flows: A proof of Ehrcnfcst time. Asymptotic Analysis 21:149-160 
(1999). 

C. Bardos, J. Garnier & G. Papanicolaou. Identification of Green's Functions 
Singularities. Inverse Problems. 24:015011 (2008). 

A. Bouzouina & D. Robert. Uniform semiclassical estimates for the propagation 
of quantum observables. Duke Math. Jour. lll(2):223-252 (2002). 

F. Brcnguier, N. M. Shapiro, M. Campillo, V. Fcrrazzini, Z. Duputel, O. 
Coutant & A. Ncrcessian. Towards Forecasting Volcanic Eruptions using Seis- 
mic Noise. Nature Geoscience doi:10.1038/ngeol04 (2008). 

M. Campillo & A. Paul. Long range correlations of seismic codas. Science 
299:547-549 (2003). 

Y. Colin de Verdiere. Methodes semi-classiques et theorie spectrale. Lecture 



Notes in Preparation http://www-fourier.ujf-grenoble.fr/~ycolver/ 



Y. Colin dc Verdiere. A semiclassical inverse problem II: reconstruction of the 
potential. ArXiv: |math-ph/0802164ft . 

A. Derode, E. Larosc, M. Tantcr, J. dc Rosny, A. Tourin, M. Campillo & M. 
Fink. Recovering the Green's function from field-field correlations in an open 
scattering medium. J. Acoust. Soc. Am. 113:2973-2976 (2004). 

A. Derode. E. Larose, M. Campillo & M. Fink. How to estimate the Green's 
function of a heterogeneous medium between to passive sensors? Application 
to acoustic waves. Applied Physics Lett. 83(15):3054-3056 (2003). 

M. Dimassi & J. Sjostrand. Spectral Asymptotics in The Semi- Classical Limit. 
London Math. Soc. Lecture Notes Ser. 268 (1999). 

J. Duistermaat. Fourier Integral Operators. Birkhduser, Boston, 1996. 

T. Duvall, S. Jefferies, J. Harvey & M. Pomerantz. Time-distance in helioseis- 
mology. Nature 362:430-432 (1993). 

C. Emmrich & A. Weinstcin. Geometry of the Transport Equation in Multi- 
component WKB Approximation. Commun. Math. Phys. 176:701-711 (1996). 

L. Evans & M. Zworski. Lectures on semiclassical analysis (Version 0.2). Avail- 



able at http://math.berkeley.edu/~zworski/ 



I.M. Gel'fand & NY. Vilcnkin. Generalised functions, vol. 4: applications to 
harmonic analysis. Academic Press, 1964. 

P. Gouedard, L. Stehly, F. Brenguier, M. Campillo, Y. Colin de Verdiere, E. 
Larose, L. Margerin, P. Roux, F.J. Sanchcz-Sesma, N. Shapiro & R. Weaver. 
Cross-correlation of random fields: mathematical approach and applications. 
Geophysical prospecting 56:375-393 (2008). 

L. Hormander. The spectral function of an elliptic operator. Acta Mathematica 
12:193-218 (1968). 

O.I. Lobkis & R.L. Weaver. On the emergence of Green's function in the 
correlations of a diffuse field. J. Acoust. Soc. Am. 110:3011-3017 (2001). 



32 



[19] J. Rauch & M. Taylor. Decay of Solutions of Non-dissipativc Hyperbolic Sys- 
tems on Compact Manifolds. Comm. Pure Appl. Math. 28:501-523 (1975). 

[20] M. Reed & B. Simon, Methods of modern Math. Phys. III. 

[211 Ph. Roux, K. Sabra & W. Kuppcrman. Ambient noise cross correlation in free 
space: Theoretical approach. J. Acoust. Soc. Am. 117:79-84 (2005). 

!l K. Sabra, P. Roux, A. Thodc, G. D'Spain, W. Hogliss & W. Kuperman. Using 
Ocean Ambient Noise for Array Self-Localization and Self-Synchronization. 
IEEE J. of Oceanic engineering 30:338-346 (2005). 

[231 F. Sanchez- Sesma, J. Perez-Ruiz, M. Campillo & F. Luzon. The elastodynamic 
2D Green's function retrieval from cross-correlation: the canonical inclusion 
problem. Geophysical Research letters 33:13305 (2006). 

[24] N. Shapiro & M. Campillo. Emergence of broadband Rayleigh waves from 
correlations of the ambient seismic noise. Geophys. Res. Lett. 31:L07614 (2004). 

[25] N. Shapiro, M. Campillo, L. Stehly & M. Ritzwoller. High Resolution Surface 
Wave Tomography From Ambient Seismic Noise. Science 307:1615 (2005). 

[261 R- Schubert. Scmiclassical behaviour of expectation values in time evolved 
Lagrangian states for large times. Commun. Math. Phys. 256:239:254 (2005). 

[271 L. Schwartz. Radon measures on arbitrary topological spaces and cylindrical 
measures. Oxford University Press, 1973. 

[281 J- Sjostrand. Projecteurs adiabatiques du point de vue pseudo-diffcrentiel. C. 
R. Acad. Sci. Paris, Serie 1, 317:217-220 (1993). 

[291 S. Teufel. Adiabatic perturbation theory in quantum dynamics. Lecture notes 
in mathematics 1821 (2003). 

[301 F. Treves. Introduction to pseudo differential and Fourier integral operators. 
Plenum Press, New York, 1980. 

[311 R.L. Weaver & O.I. Lobkis. On the emergence of Green's function in the 
correlations of a diffuse field: pulse echo using thermal phonon's. Ultrasonics 
40:435-439 (2002). 

[321 R.L. Weaver & O.I. Lobkis. Ultrasonics without a source: thermal fluctuation 
correlations at MHz frequencies. Phys. Rev. Lett. 87(13):134301-134304 (2001). 

[33] R.L. Weaver & O.I. Lobkis. Diffuse fields in open systems and the emergence 
of the Green's functions (L). J. Acoust. Soc. Am. 116(5):l-4 (2004). 

[34] R.L. Weaver. Information from seismic noise. Science 307:1568-1569 (2005). 

[35] N. Wiener. Extrapolation, Interpolation and Smoothing of Stationary Time 
Series. MIT and John Wiley, 1950. 

[36] K. Yosida. Functional Analysis. Springer- Verlag, Grundlheren der Math. Wiss. 
123, 1965. 



33 



